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Abstract 

We have used a self-consistent Skyrme-Hartree-Fock plus Continuum- RPA model 
to study the low- multipole response of stable and neutron/proton-rich Ni and Sn 
isotopes. We focus on the momentum-transfer dependence of the strength distri- 
bution, as it provides information on the structure of excited nuclear states and in 
particular on the variations of the transition form factor (TFF) with the energy. 
Our results show, among other things, that the TFF may show significant energy 
dependence in the region of the isoscalar giant monopole resonance and that the 
TFF corresponding to the threshold strength in the case of neutron-rich nuclei 
is different compared to the one corresponding to the respective giant resonance. 
Perspectives are given for more detailed future investigations. 

PACS: 21.60.Jz; 21.10.Re; 24.30. Cz 

1 Introduction 

The multipole response of nuclei at energies below the quasielastic regime is character- 
ized by various collective excitations. Among these, giant resonances (GRs) of stable 
nuclei have been investigated for decades, both theoretically and experimentally ^ElEl- 
Macroscopically, they are interpreted as small-amplitude excitations corresponding to 
the propagation of zero sound in the nuclear medium. Collective modes of excitation at 
energies below GRs include surface vibrations and rotations of the nuclear body. Trans- 
verse modes of excitation, corresponding to oscillations of the current distribution, also 
exist - a well-known example being the twist mode '3' and perhaps the toroidal mode 
[S]. At energies higher than GRs, overtones of GRs are examined [3 El EI- Despite the 
long-time effort devoted to the understanding of nuclear collective motion, several issues 
remain open; as far as the electric response of spherical stable nuclei is concerned, we 
may mention as examples the fine structure of GRs, the experimental energy-weighted 
sum rule (EWSR) of compression modes (isoscalar monopole and dipole resonances), 
the relation between the energy of the latter and the value of nuclear matter incom- 
pressibility, the identification of isovector quadrupole GRs and the nature of low-lying 
collective modes 1312]. 

Nowadays, it is feasible in addition to study experimentally the properties of unsta- 
ble nuclei, in radioactive ion beam facilities. Access to their multipole response appears 
possible using inverse kinematics or double storage rings. The physics of nuclei far 
from stability is rich in new phenomena, owing to the presence of very weakly bound 
nucleons, the unusual ratio N/Z (for given mass number) and the large difference be- 
tween the proton and neutron Fermi energies. The structure of the particle continuum 
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and the nature of low-lying strength in unstable nuclei is of particular importance, as 
they provide information on decay properties, polarizabilities, shell structure etc. At 
present, limited experimental information is available. The giant dipole resonance of 
neutron-rich Oxygen isotopes has been under study, as well as the dipole response of 
various light unstable nuclei (eg. neutron-rich He, Li, Be, C isotopes and proton-rich 
®Be and ^'^O )tlO|. Information is available on the low-lying quadrupole strength of 
various unstable nuclei Theory has been anticipating experiment to some extent, 
and has by now provided interesting predictions regarding the low-energy response and 
GRs of nuclei far from stability. For example, a considerable mixing of isoscalar and 
isovector transitions and an increased fragmentation of GRs are expected. A special 
feature of the response of very neutron-rich nuclei is the so-called threshold strength, 
namely the significant amount of isoscalar and non-collective strength predicted to oc- 
cur just above the neutron threshold. In the case of open-shell nuclei, the shell model 
[E] and the Quasiparticle RPA (QRPA) [131 El CSl HEI have been apphed. Relativistic 
RPA (RRPA) has been used as well ^1 El ■ Results on doubly magic nuclei have been 
obtained with standard RPA (with or without inclusion of higher-order configurations) 
or Continuum RPA (CRPA) methods; in particular, the self-consistent Skyrme-Hartree- 
Fock -t- Continuum RPA (SHF-I-CRPA) has been used for a systematic study of the 
response of nuclei far from stability by the group of I. Hamamoto ^1 OOl 1^ 1211 1221 
and recently by our group |24II25| . For more information on the relevant theoretical and 
experimental activity, one can refer, for example, to the contributions in Ref. 

The RPA method describes satisfactorily the transition densities of collective states 
in stable nuclei. In the past the RPA has also been used to study transition densities 
and currents of individual excitations in the case of unstable nuclei and to compare 
with the behaviour of stable nuclei [TniEniiniE21l2SlElESlEElEIlEHl- In this 
work we focus on the transition densities and form factors corresponding to electric 
excitations of stable and unstable nuclei, adopting an approach which allows a systematic 
study. Using the SHF-f CRPA method, we will examine how the transition strength 
distribution of selected Sn and Ni isotopes, and for low multipolarities, varies with the 
momentum q transferred to the system. For this we will consider an external field of 
the form jLili")- The momentum dependence of the nuclear response can be studied 
experimentally using inelastic electron scattering. Further information accessible in 
inelastic electron scattering experiments concerns the transition current density, which 
has been the subject of a separate work |25| . 

The purpose of the present pursuit is twofold. First of all, it is important to have a 
reliable microscopic description of the transition density characterizing a nuclear excita- 
tion - either of a stable or an unstable nucleus - not only to gain insight into the nature 
of the excitation, but also because theoretical transition densities enter the analysis of 
electron- or hadron-scattering experiments aiming at identifying GRs. Theoretically, the 
transition density associated with a particular type of excitation approaches hydrody- 
namical behaviour in the case of very collective states, such as GRs, but in general it 
is expected to be energy dependent. Therefore, macroscopic descriptions may be inade- 
quate in certain cases. Our approach will help us to pinpoint cases where the transition 
density varies significantly with excitation energy. Possible candidates are resonances 
with large fragmentation width. Second, we wish to contribute to the theoretical under- 
standing, currently being formed, of the response of unstable nuclei, as a guidance to 
future experimental activity. Inevitably, some of the results obtained using conventional 
models and interaction parametrizations (fitted to describe the properties of stable nu- 
clei) will be proven inaccurate by the future experimental research, but there is much 
to be learned in the process. 

In Sec. 2 the quantities of interest are defined and our method of calculation is 
presented. Results for spherical, closed-shell Sn and Ni isotopes are presented and 
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discussed in Sec. 3. Conclusions and perspectives are given in Sec. 4. 



2 Definitions and method of calculation 

We consider the response of spherical, closed- (sub)shell nuclei to an external field of the 
type 

A 

i=l 

where the variable Ti = p or n labels the isospin character -proton or neutron- of the z-th 
particle. For an isoscalar (IS) field, VL{r,p) = VL{r, n) = VL{r) and for an isovector (IV) 
one VL{r,p) = —VL{r,n) = VL(r). For L = 1 we use an effective charge equal to — ^ 
for protons and for neutrons. In the following, the isospin label will be suppressed for 
the sake of simplicity. 
We set 

VL{r) = [47r(2L+l)]i/2j^(gr), 

where Jl is a spherical Bessel function. In the long-wavelength limit qR — > (i? is the 
nuclear radius), we obtain the usual multipole operator of the form , where K — L 
for L > and X = 2 for L = 0. The transition density 5plo{t,E), characterizing the 
excited natural-parity state \LQ) of energy i?, is determined by its radial component 
SpL{r, E), where 

5pUr. E) - {2L + l)-^'^5pL{r, E)YLo{e, c^)/r\ (1) 

For IS (IV) transitions, this is the sum (difference) of the proton- and neutron-transition 
densities. The strength distribution 

S{E,q)^Y.\^mf)?6{E-Ef) 
f 

(where |/) are the final states, excited by the g— dependent external field V and Ef 
their excitation energies) is related to the transition form factor FL{q'^,E), i.e. to the 
Fourier transform of SpLQ{r, E). In particular, since we are dealing with continuous 
distributions, we write the strength in a small energy interval of width AE at energy E 
as 

S{E,q) = ^""^^^^ I r dr6pL{r.E)3dqr)\^ 



(2L + 1) 
AE 



dVVio(r,i?)e'«-T « \FL{q\E)\\ (2) 



As long as charge-current conservation is respected, the calculated S{E, q) may be in- 
terpreted as the longitudinal response function entering the analysis of inelastic electron 
scattering within the Born approximation j29| . 

Our aim is to search for variations of the transition density in the region of collective 
excitations. By looking at the strength distribution S{E, q) corresponding to different 
values of momentum q, we will "probe" the various momentum- (Fourier-) components 
of SpL, according to Eq. (2). In a scenario where the SpL (and Fl) turns out to be 
the same, within a proportionality factor f{E), for all transitions making up a strength 
distribution S{E,q), the shape of the calculated S{E,q), determined by f^{E), should 
not change with q. In realistic cases there will be non-trivial variations of the SpL with 
the energy, revealing the different microscopic structure of the various transitions. 
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Following the standard RPA method, we consider particle-hole (ph) excitations, built 
on top of the Hartrec-Fock (HF) nuclear ground state and subjected to the ph residual 
interaction (HF+RPA method). In particular, the quantities introduced above are cal- 
culated using a Skyrme - Hartree-Fock (SHF) plus Continuum - RPA (CRPA) model. 
For the HF ground-state, the numerical code of P.-G. Reinhard 'SH' is used. The cal- 
culation of the response function (unperturbed HF, as well as RPA) is formulated in 
coordinate space, as described in 31 , 32 , 33 , 34 , 35 . The radial part of the unperturbed 
ph Green function, of multipolarity L and specified isospin character, is given by: 

C'l^^.^/.^y_y ( ip\\OL\\h};{p\\0'L\\h),.' ^ {h\\0'L\\p):,{h\\OL\\p)r \ ^ 

ph ^ 

where Ol (or O'^) is one of the operators Yl, [11 ® (V^ -h V'^)]l, [Fiii ® (V - V')]l and 
[Yl±i ® (V -I- V')]l. The sign of the second term depends on the symmetry properties 
of the operators Ol and under parity and time-reversal transformations. Spin oper- 
ators have been omitted in the present calculation. With h (p) we denote the quantum 
numbers of the HF hole (particle) state and Eph = £p — £h is the energy of the unper- 
turbed ph excitation. The particle continuum is fully taken into account, as described 
in |31l 1341 13()j . A small but finite ImE ensures that bound transitions acquire a finite 
width 

The RPA ph Green function is given by the equation 

Gr^ = [l + GlV,,r'Gl, (4) 

which is solved as a matrix equation in coordinate space, isospin character and operators 
Ol- The ph residual interaction Kos is zero-range, of the Skyrme type, derived self- 
consistently from the Skyrme-HF energy functional |^ E2] . From the Green function 
(^RPA _ Qi^ _ the strength distribution S{E, q) is obtained as 

S{E,q) = 4(2L + l)Im j jL{qr)Gf"'{r,r';E)jL{qr')drdr' . (5) 

Within the RPA, giant resonances (GRs) are generated as coherent ph excitations 
between major ni shells; for example, AN=2 for the giant monopole and quadrupole 
resonances (GMR and GQR) and AN=1 for the isovector (IV) giant dipole resonance 
(GDR), where N = 2n + £ is the energy quantum number of a shell. Therefore, a GR lies 
energetically in the neighbourhood of the ANfko region, where the energy quantum hco ~ 
41A~^/^ MeV, within the simple harmonic-oscillator model. Its precise position relative 
to this value (lower or higher) is determined by the ph residual interaction (attractive 
or repulsive). "Overtones" of GRs are generated as A A''' = AA^ -I- 2 excitations. An 
example is the isoscalar (IS) GDR, which lies energetically in the '^hio region. The 
HF-I-RPA method can also describe low- lying transitions with AA^=0 (2+, 4+ etc), 
provided that the nucleus considered is not n£— closed. 

3 Results and discussion 

We have applied our method to six cases of spherical, closed-shell nuclei lying inside, 
close to or far away from the valley of stability. In particular, we have examined the 
isotopes 56,78,ii0j,jj and ioo,i2oa32gj^ ^ be^-^ jg ^he next heavier Z = N, doubly closed 
nucleus after '^"Ca . It is /3— unstable, but lies very close to the stability fine; being 
closed-shell, ^^Ni is the most stable Ni isotope where our model could be apphed, since 
pairing is not included in the model. ^*Ni is a neutron-rich isotope, possibly doubly 
magic As was done before in Refs. JHIEHI; we have used the extremely neutron-rich 
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configuration ^^''Ni as an academic example of a closed-sliell isotope in tfie vicinity of 
the neutron drip line. According to SHF results, the N = 82 closure may still be valid in 
the Ni region although not conclusively. Lying on the proton-rich side of the nuclear 
chart and with a half life of the order of Is [211011], ^""Sn may be the heaviest N ^ Z 
nucleus inside the proton drip line. '^^"Sn is a stable, doubly magic Sn isotope, while 
^^^Sn is a neutron-rich isotope with a half-life of 39.7s 0J and among the most magic 
heavy nuclei, as no excited states of this nucleus have been detected below 4 MeV 

We have calculated the IS and IV monopole (ISM and IVM), IV dipole (IVD) and IS 
and IV quadrupole (ISQ and IVQ) response of the nuclei 56,78, iiojs^j iio,i20,i32g^ 
ing the RPA method described in the previous section, for q = 0.2, 0.4, 0.6, 0.8, 1.0 fm~^ . 
We have employed the Skyrme parametrization SkM*'43', tailored to describe GRs of 
stable nuclei and used in previous studies of the response of exotic nuclei as well, eg. in 
Refs. [HEOIEIlEllESllHlESlEiEZI- We have also used the Skyrme force MSk7 gH, 
whose parameters were determined by fitting the values of nuclear masses, calculated 
using the HF+BCS method, to the measured ones, for 1888 nuclei with various values 
of \N — Z\l A. The two forces have similar nuclear-matter properties, except for the 
effective mass m*. The results obtained with the two Skyrme parametrizations agree 
qualitatively. Therefore, we only show results obtained with SkM*. 

Selected results are presented in Figs. lllSl In Figs. I6I8I the radial part of the proton- 
and neutron- transition density is plotted for some cases. We notice that, in all examined 
cases, as q increases, the strength distribution is shifted to higher energies, which can be 
interpreted as the onset of the quasielastic peak. The continuum becomes increasingly 
important. Also, overtones of giant resonances become visible. For instance, in the ISM 
response. Fig. ^ strength is shifted from the 2?ia; region (the IS GMR) to the Ahuj 
region. In the IVD response, Fig. 01 strength is found in the ihuj region as q increases. 
Excitations of single-particle character, with density oscillations taking place in the 
interior of the nucleus, give rise to this behaviour of the form factor - cf., for example. 
Fig. |S1 right panel, where the IVM transition density of ^■^^Sn at 46 MeV is shown. In 
the medium-heavy nucleus ^^Ni the shift takes place more slowly as a function of q than 
in heavier nuclei, such as ^^"^Sn , because a narrower density distribution corresponds to 
a broader form factor. 

Next we discuss our results in more detail. 

3.1 Isoscalar monopole response and compression modes 

In Fig.^the ISM response of ^^Ni , ^-'^^Ni and ^^"Sn is presented. Observing the values 
of the strength distribution of ^^°Sn for the various values of g, we notice that the form 
factor of the weak peaks on the high-energy tail of the GMR has a weaker g— dependence, 
between q = 0.4 — 0.8 fm"-'^, compared to the GMR peak. The same holds for the other 
Sn isotopes and for ^^Ni (not shown). The GMR width of the medium- heavy nucleus 
^^Ni is large, compared to the GMR width of heavier nuclei. In Fig. ^ (left panel) two 
distinct energy regions can be recognized in the ISM strength distribution of ^^Ni : the 
region P< below £q « 23 MeV, and the region P> above £q. According to macroscopic 
models, the GMR is a uniform compression mode whose transition density has a node 
at the nuclear surface. The node would then occur at radius R ~ \.2A^/^ = 4.6 fm 
(or 4.3 fm, if we use our SHF result for the radius) in ^^Ni . Therefore, the transition 
density would show maximal overlap with the function jo{qr) (whose 1st root equals tt) 
for q = qoi ~ n/R — 0.68 (or 0.73) fm"-'^. It seems that the form factor in the region 
P< follows, at least approximately, this type of behaviour, since it reaches a maximum 
between 0.6 and 0.8 fm"-'^. The form factor in P> is maximized at a larger value of q 
and therefore it does not correspond to such a picture. Indeed, as we observe in Fig. El 
the transition density in the regime P> has a node at a smaller radius than in P<, by 
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about 0.5 fm, and one more node at a larger radius. 

The fragmentation of the GMR is a typical feature of light and medium-heavy nuclei. 
It seems therefore important, for an accurate evaluation of the EWSR and the centroid 
energy of the GMR in those nuclei, to take into account the energy dependence of the 
transition density in the analysis of the relevant a— scattering experiments - see also the 
analysis in Ref. 021 • The same may hold for the IS giant dipole resonance (GDR), which, 
in addition, appears to have large width even for nuclei as heavy as ^°*Pb . Moreover, 
the amount and nature of the strength detected below the main IS GDR peak (which 
is located at approximately 23 MeV for heavy nuclei) has not been clarified. For a 
recent report on the issue see eg. the contribution of U. Garg in Ref. Given that 
the properties of the IS GMR and GDR are used for determining the value of nuclear 
matter incompressibility, a detailed examination may be recommended. We have not 
presented calculations of the IS GDR here, as our method lacks full self-consistency (due 
to the omission of spin-depended terms from the residual interaction) and therefore our 
results would not be free of spurious components. 

3.2 Isoscalar quadrupole response 

The two peaks in the ISQ strength distribution of the nuclei considered here are the low- 
lying collective 2+ transition (first peak) and the IS GQR (second peak). In the cases 
of '^^Ni (see Fig. 12), ''^Ni (not shown) and ioo,i20gj^ (^-^q^- gi^Q-^^i^)^ the two peaks show 
similar behaviour as a function of g up to (7 « 0.8 fm^^. This is due to the fact that in 
these nuclei, the transition density of both states is peaked close to the surface. It should 
be noted, however, that the nature of the two peaks may be quite different, as the low- 
lying state is expected to be characterized by non-neglig ible vorticity [21EH1- We also 
find that the transition density corresponding to the low-lying state has a node inside the 
nucleus. In the case of ^'^■^Sn (Fig. the form factor of the first peak does not follow 
the behaviour of the GQR form factor. In ^^"Ni (Fig.|2l, there is a clear difference. The 
low-energy peak loses its strength faster than the GQR, as q increases, behaving like 
the threshold strength of other multipolarities - which we discuss in § 13.41 Respective 
transition densities are plotted in Fig. [7| In ^^"Ni , a third peak occurs between the 
low-lying state and the GQR, also showing different behaviour. Its strength peaks at 
higher q than the GQR and the low-lying-state strength. It is a non-collective state with 
a transition density localized in the interior of the nucleus. Such a peak occurs also in 
the case of ^^Ni (not shown) and ^^^Sn . Corresponding transition densities are also 
plotted in Fig.[7| They indicate that these secondary peaks are not of pure IS character. 

The experimental width and fragmentation of the IS GQR can be much larger than 
accounted for by first-order RPA calculations. In order to reproduce the experimental 
width of the GQR one has to take into account higher-order configurations than Iplh 
|4()l I47j . Then one would be in a position to examine the energy dependence of the 
transition density in the region of the GQR. 

3.3 Isovector strength distributions 

As shown in Fig.|3| the IVM response of neutron-rich nuclei is dominated by two struc- 
tures, namely the broad peak of the IV GMR above 25 MeV, higher than the IS GMR 
due to the repulsive IV residual interaction, and an amount of strength at the IS GMR 
region. The latter is of IS character (see Fig.|Hl left panel) and exemplifies the admixture 
between IS and IV transitions in neutron-rich nuclei. A similar situation is observed in 
the IVQ response (Fig. |3J); most of the IVQ strength, making up the IV GQR, fies 
above 20 MeV, but an amount of strength remains in the IS GQR region. In all nuclei, 
the IV GMR and GQR have a width of several MeV. 
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No significant energy dependence of the transition density is observed in the region of 
the IV GMR and GQR for the low values of q examined here. With the exception of the 
threshold strength in very-neutron-rich nuclei (see ij l3.4|l . the IVM strength distributions 
in Fig. 01 for q = 0.2,0.4 and 0.6 fm~^ are similar to each other. The same seems to 
hold for the IVQ distributions (Fig.EJ, in the region of the IV GQR, in spite of the rich 
fine structure of the latter. In the case of the IVD distributions in Fig. ^ this holds to 
a lesser extent. 

The nature and systematics of the low-lying dipole strength is of particular impor- 
tance in connection to astrophysical processes. A detailed examination of the IV as well 
as IS dipole strength distribution and also the corresponding vorticity strength distri- 
bution for various nuclei, including the momentum dependence of these distributions, 
using a fully self-consistent method, should be able to clarify - from a theoretical point 
of view - important issues such as the admixture of IS and IV transitions at low energies, 
the existence of toroidal modes and the energy dependence of the transition densities 
and currents in the region of the IS GDR. 

3.4 Threshold strength 

As mentioned in the introduction, the response of very neutron-rich nuclei is character- 
ized by the so-called threshold strength [191 1201 1^ I23L I27| . The effect is expected when 
the neutron threshold lies much lower than the GR region - which begins typically at 
about 10 MeV. The values of the particle threshold energy Et for the isotopes 56,78, iiO]\jj 
and ioo,i20,i32gj^ ^^-^^^ ^j^g multipolarilties considered here are presented in Fig. 
According to our results, the threshold strength vanishes as q increases. This is seen 
clearly in the case of ^"'^^Ni . As we observe in Fig. ^ (middle panel), the ISM threshold 
strength, appearing as a "shoulder" at low energy, behaves differently than the peak at 
around 13.5 MeV; thus we have a way to discriminate between the "threshold strength" 
and the main peak - which represents the GMR. The same holds for the IVM response 
of "'^^^Sn (see Fig. |3| right panel) and similarly for the IVD (see Fig. ^ middle panel). 
The distribution of loosely bound neutrons at large distances from the nuclear center 
results in form factors Fl with maxima at low values of giving rise to this type of 
behaviour of the threshold strength. The effect of the excess neutrons and the origin 
of the threshold strength is demonstrated in Figs. |51 (^^'^Ni , compare right-top panel 
with the other panels) and |H1 (^'^^Sn , left panel), where typical examples of transition 
densities corresponding to the threshold region are presented. 

In Ref. 22 the ISM response of the neutron- rich nucleus ^°Ca is examined - where a 
significant amount of threshold strength appears as well. It was found that in the region 
of the GMR the transition density compares rather well with the Tassie model prediction, 
whereas in the threshold region it is extended at large distances and it originates mostly 
from neutron excitations. A similar effect was predicted also in the case of the ISQ 
response of ^^0 Our results are in concordance with these findings. 

4 Conclusion and perspectives 

In this work we have used the SHF-fCRPA method to study the low-multipolarity 
response of selected Ni and Sn isotopes and to examine in particular variations of the 
transition density and form factor in the region of collective excitations. According 
to our results, the transition density may show considerable energy dependence in the 
region of the IS GMR. This should be taken into account in the analysis of relevant 
hadron scattering experiments. The form factor, corresponding to the threshold strength 
in very-neutron-rich nuclei, is narrower compared to the form factor of the respective 
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giant resonance. This result, owing to the losely bound neutrons outside the core, is 
independent of L. In the region of IV GMR and GQR resonances, no significant energy 
dependence is observed for q values below 1.0 fm"^. A detailed examination of the 
IS and IV dipole strength distribution, as well as the corresponding vorticity strength 
distribution, should be the subject of future work. Special care should be taken of the 
spurious center-of-mass motion. Important issues to be addressed include the admixture 
of IS and IV transitions at low energies, the existence of toroidal modes and the energy 
dependence of the transition densities and currents in the region of the IS GDR. 
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Figure 1: ISM strength distribution as a function of energy and momentum transfer. 
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Figure 2: ISQ strength distribution as a function of energy and momentum transfer. 
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Figure 3: IVM strength distribution as a function of energy and momentuni transfer. 
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Figure 4: IVD strength distribution as a function of energy and momentum transfer. 
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Figure 5: IVQ strength distribution as a function of energy and momentum transfer. 
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Figure 6: The transition density (in arbitrary units) for protons (fuU hues) and neutrons 
(dashed hues), corresponding to ISM transitions of the nuclei ^^Ni (left panel) and ^^"^Ni 
(right panel) at the indicated values of the energy. 
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Figure 7: The transition density (in arbitrary units) for protons (full lines) and neutrons 
(dashed lines) corresponding to the first (collective) and second low-lying ISQ peaks and 
to the IS GQR, for the nuclei ^«Ni , "^Ni and ^^'^Sn . 
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Figure 8: The transition density (in arbitrary units) for protons (full linos) and neutrons 
(dashed lines) corresponding to IVM transitions of the nucleus ^''^Sn at the indicated 
values of the energy. 
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Figure 9: The calculated particle threshold energy Et plotted vs the mass number A, for 
the isotopes 56,78,ii0]sjj ^^^^ ioo,i20,i32gjj g^jjjj £qj. i = 0, 1, 2. Lines connecting isotopes 

of the same element are drawn to guide the eye. 
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